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Abstract 

Centre manifold techniques are used to derive rationally a description of the dynamics 
of thin films of fiuid. The derived model is based on the free-surface ri(x,t) and the 
vertically averaged horizontal velocity u(x,t). The approach appears to converge well 
and has significant differences from conventional depth-averaged models. 

1 Introduction 

Consider the dynamics of an incompressible viscous fiuid fiowing down an inclined plane 
of angle 6. Long wave, low Reynolds number dynamics of the fiuid film may be described 
by a well known equation [3, eq.(l)]. In the limit of small departures from a fiat in- 
terface, it reduces to a Kuramoto-Sivashinsky equation [3, eq.(2)] or [4, eq.(20)]. Such 
long-wave models may be put on a firm theoretical footing with the aid of centre mani- 
fold techniques as has been done for Taylor dispersion[9, 19], fiow reactors[l] and beam 
theory [17]. However, in many applications these particular long- wave models of fiuid 
films have limited usefulness and a two-mode long-wave model of the dynamics, ex- 
pressed in terms of the film thickness rj and mean velocity u, is preferred [3, pllO]. This 
note shows how to use centre manifold techniques to derive such two-mode long-wave 
models. 

An advantage in using centre manifold techniques is that in a range of applications 
there is rigorous support for the resulting low-dimensional models [2]. However, rigor- 
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ous theory does not yet cover this application. Nonetheless, based upon a large spectra 
gap. Section 2, I argue that a useful two-mode model can be constructed. By modify 
ing the Navier-Stokes equations. Section 3, I fit the problem into the centre manifoL 
formalism and seek an asymptotic solution. The leading-order model of the dynamics 
Section 4, exhibits significant quantitative differences with other models based upoi 
depth-averaging (or equivalently depth integration); these differences persist in a mor 
careful examination of the convergence of the model. Section 5. Proceeding to second 
order in the analysis. Section 6, I derive a dynamical model which includes terms pro 
portional to the Reynolds number; these terms are entirely absent in depth-average( 
equations. Consequently, I deduce that depth-averaging is qualitatively and quanti 
tatively unsound except perhaps for low Reynolds number fiows. Instead, I propos' 
this centre manifold approach for developing models that agree quantitatively with th 
original governing differential equations. 

With the a;-axis along and the y-axis perpendicular to the plane of a bed at angle 6 
the Navier-Stokes and continuity equations for the dynamics of a film of incompressibl 
fiuid are: 
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3 THE CENTRE MANIFOLD ANALYSIS 
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where I have scaled quantities so that fluid density p = 1, and typical magnitudes of 
fluid thickness and horizontal velocity are of size 1. Thus, in terms of a typical velocity 
U and a typical fllm thickness H , the dimensionless constants appearing in the above 
equations are as follows: 

• g is the reciprocal of the Froude number l/Fr = gH/U^; 

• i^ is the reciprocal of the Reynolds number 1/ Re = jj,/ HU p; 

• and (T is the Weber number We = T/pU^H . 

There is a one-parameter family of equilibria of the governing equations (1-7), they 
correspond to uniform shear flow. In terms of a constant h an equilibria is rj = h, 
p = —g cos6(h — y)^ u= ^ '^™ {hy — |j/^), and v = Q. Using centre manifold techniques 
we can construct low-dimensional dynamical models based on these equilibria. Note 
that because there is a family of equilibria, essentially parameterised by the free-surface 
height rj, the derived model will permit order 1 changes in rj over the flow domain — just 
provided that the surface variations are slow. 

2 The spectrum 

In the absence of any x variation, the spectrum of perturbations to any one of these 
equilibria is {Aj} where: 

• Ao = corresponds to the freedom to vary the fluid depth h; 
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for n = 1 , 2 



, ^, . . . corresponding to viscously dissipating modes 



of horizontal shear m„ = sin "^ " — -y\ 



The eigenvalue permits the construction of a centre manifold model of flows which vary 
slowly in space[13]. Since, in the absence of spatial variations, there is one eigenvalue, 
then the model would be written in terms of one slowly-varying mode. The evolution 



equation of such a model would be the traditional long-wave equation [3, eq.(l)] for th 
fluid depth rj (possibly with higher-order modiflcations). 

However, there is also a clear spectral gap between Aq = allied with Ai = —i^ {jj,) 

and the other eigenvalues headed by A2 = —9^ {ji,) ■ Hence a two-mode model base( 
on the dynamics of rj and the gravest horizontal shear mode ui should be of interest ii 
many applications [15]. The model should resolve accurately dynamics and transient 
on time-scales slower that I/A2. 



3 The centre manifold analysis 



I make the gravest mode, ui 
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This modiflcation to the right-hand side shifts the eigenvalues {Ai, A2, . . .} by just th 
coefficient of the new term in u. Treating 7 as small in the asymptotic scheme (fo 
example, by the standard trick of appending dj/dt = to the set of equations, as i 
done to unfold bifurcations [2]), ui(y) corresponds to a eigenvalue, and we construe 
a two-mode centre manifold model. Ultimately, I set 7 = 1 to recover approximat 
expressions for the original problem. The model describes the long-term evolution c 
thin fllm dynamics in terms of the two critical modes rj and uui(y), slowly varying in ; 
and t. 

Having identifled the critical modes, the subsequent analysis is straightforward a 
explained in more detail for other systems elsewhere[13, 9, 17, for example]. Collectin, 
the unknown flelds into u = (u,v,p) we seek a low-dimensional centre manifold givei 
by the asymptotic expansion 



u(x,y,t) = Y(r],u,j,y) 
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where all the x and t dependence is solely through the evolution of ri(x,t) and u(x,t] 
Thus we pose an evolution of 
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4 TEE LEADING ORDER MODEL 
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In these asymptotic expansions the superscript "*'" denotes a term which is of order m 
in 7, as is explicitly shown, and implicitly of order n in the bed-slope 6 and spatial 
derivatives d/dx.[lS\ Further, it is convenient (although not essential) to consider the 
surface tension parameter a to be large, of order —2, so that surface tension effects are 
promoted to the leading order. 

Substituting this ansatz into the governing differential equations (2-8), which may be 
written in the abstract form /C^ = jCu + N(u) where C is the singular leading-order 
operator, we aim to solve 
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Then substituting the asymptotic expansions into this equation and equating terms of 
the same order gives a hierarchy, in order m and n, of equations to solve. The equations 
are made well-posed by requiring that the "amplitudes" of the parameterisation of the 
model are precisely the free-surface height rj and the mean horizontal velocity 

1 r 

u{x, t) = — I udy . 
^ Jo 

The resulting equations were solved and checked using the reduce computer algebra 
package. 

Note that I apply only formal centre manifold techniques; the rigorous theory for 
this situation has not yet been developed. The difficulty primarily lies in the "infinite" 
dimension of the centre manifold — "infinite" because it is parameterised by continuous 
functions of x. For example, the theory of Gallay[6] on infinite dimensional centre 
manifolds requires that the nonlinear perturbations are bounded, but here the nonlinear 
terms, N(u), involve the spatial derivative d/dx which is unbounded. I anticipate that 
in time theory will be developed to rigorous support these systematic techniques. 



4 The leading order model 

The leading order problem in the hierarchy of equations is (where superscripts 
omitted for clarity) 
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with the basic solution of a shear fiow in fiuid of depth rj: 
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Higher orders in the analysis show how space variations in rj and u couple to the physica 
processes in the Navier-Stokes equations to produce modifications of this leading orde 
structure together with time evolution of rj and u. 

For example, the next order adjustments to the pressure and velocity fields are, ii 
terms of the scaled vertical coordinate ( = 7rj//(2ry), 
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The three structure functions appearing in m*^'^ are plotted in Figure 1 as a function of ( 
Observe that they all have much the same shape that will tend to fiatten the horizonta 
velocity profile m'^''^(j/) (say positive): in thicker fiows due to the r\^ jv factor; if th 
free-surface slopes downwards due to the —r\x factors, as is appropriate for a convergin, 
fiow, or if the mean velocity is increasing due to the uu^ factor; or through surfac' 
tension effects, ex (Trjxxx, if the curvature increases in x. Indeed, this last effect ma; 
predict the recirculation regions seen in Figure 7(b,c) from the numerical simulation 
of Chang et al [4], especially as the predicted effect increases with Reynolds numbe 
{l/v). All of these dynamic modifications to the velocity profile are proportional to th 
Reynolds number. 

As you can imagine, the formulae for the fields rapidly become more complicated a 
higher orders. I do not record any of these details here, instead I concentrate upon th 
evolution equations for the model's parameters, rj and u. 

As usual, the above adjustments are found by inverting the singular linear operato 
C. This is only possible if the right-hand side is in the range of the singular operator 
and this solvability condition determines the evolution terms GjJ'^ and G^'^. Similar! 



4 TH^LEADING ORDER MODEL 
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Figure 1: profiles of tlie three structure functions appearing in u°'^, tlie first correction 
to the leading order horizontal velocity profile u°'° , as a function of the scaled vertical 

coordinate ( = 'iTy/(2'ri): , the uif'Ux/v term; , the —u^rjrix /v term; , the 

remaining term. 



at first order in 7 we determine G„' and G^' . The only non-zero G™'" at any order i 



G''l'^ , which gives the exact result 
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at all orders. This is a direct consequence of conservation of fiuid. 

Henceforth, I concentrate on the evolution equation for the mean horizontal velocit; 
u. The leading order evolution is governed by 
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upon setting 7 = 1 to recover a description relevant to the original problem. Comparin 
the leading order coefficients with those from equation (19) of Prokopiou et al [11]: 

• from the above uUx term the shape factor used in other approaches should be 5/4 
whereas from the above u^rjx/ri term it should be 7/6 — the value of 1.2 used b; 
Prokopiou et al is perhaps a reasonable compromise; 



• the coefficient of viscous decay, the vujr\^ term, is 7r^/4 
than the conventional 3; 



2.46, some 17% weake 



• but a more significant difference is that gravity and horizontal pressure gradient 
are actually less effective than depth-averaging indicates (by nearly 20%, the facto 
8/7r2 f« 0.81 instead of 1). 

The first two differences could be attributed to the different "shape" of the leadin, 
order velocity profile mi(j/), here sinusoidal instead of parabolic. However, the last i 
fundamentally new. The phenomenon is due to the response of the fiuid, primaril; 
sin(7rj//2/ry), being at an angle to the forcing 1 (either due to gravity or horizonta 
pressure gradients) when considered in the space of functions on [0,?y]. Consequentlj 
the forcing is less effective. 

A further comparison is made with the lubrication approximation, see equation (2) c 
Moriarty ti al [10] for example. The lubrication approximation arises directly from (11 
by assuming that the horizontal velocity is small and has viscous dissipation balance( 
by production from the x component of gravity and surface gradients. Thus 
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6 SECOND ORDER MODEL 



The coefficient 32/ ir^ is within 1.5% of the 1/3 needed to agree with lubrication theory. 
However, in order to obtain the leading order approximation (11), I introduced j to 
modify the original problem. It is necessary to re-examine the above comparisons in a 
more accurate solution of the original problem by proceeding to higher order in j. 

5 Convergence 

Computations indicate that the dependence upon j converges very quickly; or at least 
is an asymptotic series with excellent low-order summability. Solving for terms up to 
order 7®, I obtain the following for the evolution of u: 
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It seems reasonable to conclude from these expressions that the leading order in 7 is 
generally accurate to nearly two decimal places; whereas evaluating the above to terms 
in 7^ for 7 = 1 should give about four decimal place accuracy in the coefficients of the 
evolution equation. 

I conjecture that the reason this expansion in 7 works so well here is that the gravest 
eigenmode is well represented by the critical mode in the modified equation; in fact 
they are identical. Thus, here the approximation in 7 is primarily an approximation 
in the dynamics. The approximation that Ai Ri when compared with the neglected 
dynamics, headed by A2, is very good. Indeed, the rate at which the coefficients in the 
7 expansions are decreasing is clearly about A1/A2 = 1/9. 

However, earlier work on a simple dynamical system[12] has shown that we should 
expect high-order nonlinear terms to have poorer convergence. Generally expect that a 
term in u" will involve the Taylor expansion of 1/(A2 — «Ai7) which will not converge 
at 7 = 1 for terms of order n > 9. Thus although the low orders of the model converge 



well in 7, expect that it is only asymptotic in rj and u and that we should only use ; 
truncation of less than ninth order. In practise this is not a significant restriction a 
only the first few orders are of interest. 

An outstanding feature of the evolution equation (12) is that the quantitative dis 
crepancies with traditional depth-averaged models, such as equation (19) in [11], remai: 
almost unaltered. However, the small discrepancy with lubrication theory disappears 
setting 7=1 the shown coefficients in the first line of (12) sums to within roundoff er 
ror, 2 X 10~®, of the 7r^/12 needed for exact agreement. I deduce that depth-averagin, 
is quantitatively unsound for these sorts of shear fiows. 

6 Second order model 

Having established that there is very good convergence in the parameter 7, I exten( 
the model to higher order in spatial gradients and slope. Using reduce I find that ti 
third-order in 7 and to second-order in other quantities 
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+ 77^7^ ('^■'23lr]^ur]xxxx + IMGrjurjxrjxxx - O.lO^lrj^Uxrjxxx - 0A723r]ur]l^ 
lOOi^ 
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This is certainly a complicated evolution equation; which is one reason for not proceed 
ing to any higher order. ^ The striking qualitative differences between this equation an( 
other proposed models are twofold. 



Another reason is that the algebraic details are sufficiently complicated so that the reduce prograi 
was restricted to this order by computational resources. 



7 REMARKS 



• Most obvious are all the terms proportional to the Reynolds number (1/j^) in the 
above equation (13); they have no counterpart in the depth-averaged models (such 
as [11, eq.(19)]). Such terms arise through the subtle interplay of nonlinear inter- 
actions and velocity shear with vertical diffusion (they are somewhat akin to the 
effective diffusion term of Taylor's description[18] of shear-enhanced dispersion in a 
channel) — they are very hard to derive without the systematic approach of centre 
manifold theory. Being proportional to the Reynolds number they are likely to be 
of interest in the faster flows, though the factor of 1/100 indicates that the Reynolds 
number will need to be comparable to 100 before these terms are significant. 



• The other qualitative difference is the absence terms in (Tr/ri'^rixxx and cririxii^^ as 
appear in [11, eq.(19)]. Such terms would arise at the next order in this analysis 
as they involve five spatial derivatives and a is of order —2. The centre manifold 
approach suggests that other terms are more important. 

For a quantitative comparison of the other terms, with equation (19) in [11], I re- 
cast (13) in terms of the fiuid fiux q(x,t) = rju. 
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The earlier discussed discrepancies remain. In addition, instead of the terms 
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the centre manifold analysis recommends (to one decimal place) 

3.bqr]l ^AqxTjx -i.lqrjx 
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T] 
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Using these new coefficients should improve agreement between the one-dimensiona 
model and the full Navier-Stokes equations. 



7 Remarks 

Equation (13) (or equivalently (14)), in conjunction with the continuity equation (10 
form a consistent one-dimensional model for the dynamics of a fiuid film. Such centr^ 
manifold models are naturally of mixed order and so there is a great deal of fiexibilit; 
in choosing which terms to retain in any given application of the model — a decisioi 
that can be delayed until application, and based upon the particular parameters c 
the application, rather than forced by the restrictive scaling principles of primitiv 
mathematical methods such as that of multiples scale. Also, the geometric picture c 
centre manifold theory will lead to, in further work, expressions for the correct modellin 
of initial conditions [14], forcing [5] and boundary conditions [16]. Here I just focuse( 
on a derivation of the evolution equations. 

Due to isotropy, the generalisation of equations (13) and (14) to the dynamics of a thii 
sheet of fiuid is straightforward: the mean horizontal velocity and horizontal derivative 
just turn into the corresponding vector velocity and gradient; and the gravity factor 
become appropriate dot products. 

A very similar approach to the one taken here could be used in a wide range c 
physical problems such as the dynamics of slender jets [7], roll waves in muddy fiuid [8] 
and long waves on turbulent fiow (in preparation). 
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